Scaling of fluctuations in a colloidal glass 
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We report experimental measurements of particle dynamics in a colloidal glass in order to un- 
derstand the dynamical heterogeneities associated with the cooperative motion of the particles in 
the glassy regime. We study the local and global fluctuation of correlation and response functions 
in an aging colloidal glass. The observables display universal scaling behavior following a modified 
power-law, with a plateau dominating the less heterogeneous short-time regime and a power-law 
tail dominating the highly heterogeneous long-time regime. 
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I. INTRODUCTION 

Mostly due to the enormous practical importance of 
glassy systems there has been a vast literature describing 
different theoretical frameworks for glasses, yet without a 
common theory applicable to the diverse range of systems 
undergoing a glass transition. Increasing the volume frac- 
tion of a colloidal system slows down the Brownian dy- 
namics of its constitutive particles, implying a limiting 
density, (pg, above which the system can no longer be 
equilibrated with its bath 1]. Hence, the thermal system 
falls out of equilibrium on the time scale of the exper- 
iment and thus undergoes a glass transition Q. Even 
above (f)g the particles continue to relax, but the nature 
of the relaxation is very different to that in equilibrium. 
This phenomenon of a structural slow evolution beyond 
the glassy state is known as "aging" The system is 
no longer stationary and the relaxation time is found to 
increase with the age of the system, , as measured from 
the time of sample preparation. 

This picture applies not only to structural glasses such 
as colloids, silica and polymer melts, but also to spin- 
glasses, ferromagnetic coarsening, elastic manifolds in 
quenched disorder and jammed matter such as grains and 
emulsions 

In order to investigate the dynamical properties of 
the aging regime we consider the evolution of the cor- 
relation function C{t,tw) and integrated response to ex- 
ternal fields xi^T'tw)- These quantities are separated 
into a stationary (short time) and an aging part (long 
time): C(t,t^) = Cst{t - t^) + Ci,g{t,t^) and xit^t^) = 
Xst{t — tiu) + Xag(i,iu)), where we have included the ex- 
plicit dependence on in the aging part, and t is the 
observation time. 

We use a simple system undergoing a glass transition: 
a colloidal glass of micrometer size particles, where the in- 
teractions betweenparticles can be approximated as hard 
core potentials P, S ■ The system is index matched to 
allow the visualization of tracer particles in the micro- 
scope 0. Owing to the simplicity of the system, we are 
able to follow the trajectories of magnetic tracers embed- 
ded in the colloidal sample and use this information to 



understand the scaling of the local and global fluctua- 
tions of the correlation functions. 



II. 



EXPERIMENTAL SETUP 



Our experiments use a colloidal suspension consisting 
of a mixture of poly-(methylmethacrylate) (PMMA) ster- 
ically stabilized colloidal particles (radius Op = 1.5 fim, 
density pp = 1.19 g/cm^, polydispersity 14%) plus a 
small fraction of superparamagnetic beads (radius a„i = 
1.6 /im and density pm = 1.3 gj err? , from Dynal Biotech 
Inc.) as the tracers 

The colloidal suspension is immersed in a solution con- 
taining 76% weight fraction of cyclohexylbromide and 
24% cis-decalin which are chosen for their density and 
index of refraction matching capabilities 0- For such a 
system the glass transition occurs at 0g « 0.57 — 0.58 
U 1^ ' In our experiments we consider two samples 
at different densities and determine the glassy phase for 
the samples that display aging. The main results are 
obtained for sample A just above the glass transition 
= 0.58 ± 0.01. We also consider a denser sample 
B with (/)B = 0.60 ± 0.01, although this sample is so deep 
in the glassy phase that we are not able to study the 
slow relaxation of the system and the dependence of the 
waiting time within the time scales of our experiments. 

We use a magnetic force as the external perturbation 
to generate two-dimensional motion of the tracers Q on 
a microscope stage following a simplified design of 0. 
Video microscopy and computerized image analysis are 
used to locate the tracers in each image. We calculate 
the response and correlation functions in the x-y plane. 

The experimental set up is shown in Fig. ^ and Fig. 
We use a Zeiss microscope with a 50x objective of 
numerical aperture 0.5 and a working distance of 5mm. 
We work with a field of view 194^m x 155/im by using a 
digital camera with the resolution of 1288 pixels x 1032 
pixels. We locate the center of each bead position with 
sub-pixel accuracy, by using image analysis. For the con- 
densed samples A and B, we use the low frame rate 1/3 
frame/sec, while for the dilute sample C, we record the 
images at 1 frame/sec. The long working distance of the 
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objective is necessary to allow the pole of the magnet 
to reach a position near the sample. An example of the 
images of the tracers obtained in sample A is shown in 
Fig. ^ where we can see the black magnetic tracers em- 
bedded in the background of nearly transparent PMMA 
particles. An example of the trajectory in the x-y plane 
of a tracer diffusing without magnetic field is shown in 
Fig. nji. We note that this particular tracer moves away 
from two cages in a time of the order of 4 hours. 

The magnetic field is produced by one coil made of 
1200 turns of copper wire. We arrange the pole of the coil 
perpendicular to the vertical optical axis, and generate a 
field with no vertical component. Thus, the tracers move 
in the x-y plane with a slight vertical motion which is 
generated by a density mismatch between the tracers and 
the background PMMA particles. This vertical motion 
is very small at the high volume fraction of interest here, 
and therefore we calculate all the observables in the x-y 
plane. 

The magnetic force is calibrated for a given coil current 
by replacing the suspension with a mixture of 50 : 50 
water-glycerol solution with a few magnetic tracers. The 
distance between the top of the magnetic pole and the 
vertical optical axis is always fixed, which means that the 
magnetic force at the local plane depends only on the coil 
current. At a given current, we determine the velocity of 
the magnetic tracers at the focal point and calculate the 
magnetic force from Stokes's law, F — GirrjamV, where 77 
is the viscosity of the water-glycerol solution, am is the 
tracer radius, and v is the observed velocity of the tracer. 
The uncertainty in the obtained force comes from: (a) 
the uncertainty in the coil current which is 1%, (b) the 
beads, which are not completely monodisperse in their 
magnetic properties, causes a 10% uncertainty and 
(c) the magnetic field is slowly decaying in the field of 
view, causing a 4% uncertainty. 

In order to investigate the dynamical properties of the 
aging regime, we first consider the autocorrelation func- 
tion as the mean square displacement (MSD) averaged 
over 82 tracer particles, C{t,tw) = {Ax'^(t,tw))/2 = 
{[x{tw + At,tm) — a;(t^, t^)]^)/2, at a given observation 
time, t — tuj + At, after the sample has been aging 
for as measured from the end of the stirring process. 
Then, we measure the integrated response function (by 
adding the external magnetic force, F) given by the av- 
erage position of the tracers, x(t, t^) = {x(t^ + At, t^) — 



III. SAMPLE PREPARATION AT U = Q 

It is important to correctly determine the initial time 
of each measurement to reproduce the subsequent par- 
ticle dynamics. We initialize the system by stirring the 
sample for two hours with an air bubble inside the sam- 
ple (see Fig. ^p) to homogenize the whole system and 
break up any pre-existing crystalline regions. Then we 
place the sample on the magnetic stage and take images 




FIG. 1: (a) Picture of the experimental setup, (b) Schematic 
illustration of the setup, (c) Detail of an image of tracer in 
sample A. (d) Trajectory of a tracer in sample A showing the 
cage dynamics over 4 hours. The circle represents the size of 
the tracer particle, 3.2/im. 
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FIG. 2: Pair-distribution functions of (a) sample A and (b) 
sample B. 



to obtain the trajectories of the tracers that appear in 
the field of view. The initial times t = = are defined 
at the end of the stirring. After measuring all the tracers 
appearing in the field of view, a new stirring is applied, 
the waiting time is reset to zero, and the measurements 
are repeated for a new set of tracers. We have analyzed 
the pair-distribution function, ^(r), and two-time inten- 
sity autocorrelation function, g2{tw,tw + At), in order to 
test our rejuvenation technique. 

First, Fig. |3 plots the pair-distribution functions of 
sample A and sample B right after the stirring proce- 
dure. We calculate the pair-distribution functions by re- 
constructing the packings from 3D confocal microscopy 
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FIG. 3: Two-time intensity autocorrelation function, g2{tw + 
At,tm), of (a) sample A and (b) sample B. The solid black 
and navy curves correspond to t„ = Os and = 4000s, re- 
spectively. The dash gray and light blue curves correspond to 
the measurements after the second stirring with = Os and 
tnj = 4000s, respectively. The fact that the curves at — Os 
and = Os coincide indicates that sample has been fully 
rejuvenated to its initial state. The inset is the same plot as 
the main figure for different ranging from Os (solid black) 
to 4000s (solid navy). 



images of size 60/imx60/imxl5/im. We find that the 
samples do not show obvious crystallized region by di- 
rectly looking either at the pair-distribution function or 
the images taken from confocal microscopy. 

For these measurements we use a Leica confocal mi- 
croscope. The PMMA particles are fluorescently dyed 
so that they are ready to be observed by confocal mi- 
croscopy. We load the samples sealed in a glass cell on 
the confocal microscope stage and use a Leica HCX PL 
APO 63x, 1.40 numerical aperture, oil immersion lens for 
3D particle visualization in order to calculate the volume 
fraction and the pair-distribution function. 

Second, we analyze (72(^10 + ^t, t^) and show that the 
sample has been rejuvenated by the stirring process. We 
plot g2{tw-\-At, tw) for 3 situations: (a) after the first stir- 
ring process (t^ = 0), (b) after the sample has aged for 
a long time {t^^ = 4000s), (c) we then apply the stirring 
again and immediately plot the autocorrelation function 
for the new initial time (tu, = 0). The plot show the 
lack of correlation for the two individual measurements 

= and iu, =0, thus demonstrating that the sample 
has been rejuvenated. Technically, we record the tempo- 
ral image sequence of sample A and sample B, and study 
the aging by calculating two images' correlation defined 
as: 



g2{tw,tw + At) = 



{I{t^)I{t^ + M)) 



(1) 



where / is the average gray-scale intensity of a small 
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FIG. 4: Scaling plot of sample A for the scaled autocorrela- 
tion, tw°'C, and scaled integrated response, kBTastw'^Xy 
a function of the time ratio, At/tm'', for different waiting 
times. The black dash line is a linear fit which indicates that 
Tcff = 690-/^. The inset is a plot of sample B for autocorre- 
lation, C, and integrated response, ksTesXy ^ ^ function of 
At at = 100s. The black dash line is a linear fit which 
indicates that Teff = 16GQK. 



box (20 pixels x 20 pixels, PMMA particles' diameter 
is roughly equal to 20 pixels), and /(t^) and I{t) come 
from the images at different times tu, and t, respectively. 
We cut one image (1288 pixels x 1032 pixels) into many 
boxes, and calculate the correlation of two boxes located 
at the same position in the two images. The average 
(...) is taken over all the boxes in one image. Fig. |3| 
plots the correlation function, g2{tw,tw + At), of sample 
A and sample B, which shows aging behavior (see the 
inset). Fig. |3| shows how g2{tw,tw + Af) calculated after 
two stirring processes at tw = and 1^=0 coincide, 
indicating that the sample can be fully rejuvenated to its 
initial state by our stirring technique. 



IV. SCALING ANSATZ FOR THE GLOBAL 
CORRELATIONS AND RESPONSES. 

Insight into the understanding of the slow relaxation 
can be obtained from the study of the universal dynamic 
scaling of the observables with Based on spin-glass 
models, different scaling scenarios have been proposed 0, 
ITl| for correlation and response functions. Our analysis 
indicates that the observables can be described as 



At, K 



(2) 



where fo and /m are two universal functions and a 
and P are the aging exponents. Evidence for the validity 
of these scaling laws is provided in Fig. Qlwhere the data 
of the correlation function and the integrated response 
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FIG. 5: (a) Scaled {x\2) as a function of scaled roi for given 
ratio ^tjt^ — 10. The linear part for different tw and At are 

coincident and well described as tZ^^{xl2) — (— 0.28)tSi^^roi 
, indicated by the dashed line, (b) Probability distribution 
function(PDF) of the scaled local correlations for the ratios 



Probability distribution function(PDF) of 



the scaled local correlations for the ratios At/t^ = 60. The 
inset is a log-log plot of averaged PDF of the scaled local 
correlations for three ratios At/t^ — 1, 10, 60, the tails are 
coincident and well fitted by power law with a index -1.25. 
(d) Probability distribution function(PDF) of the scaled local 
integrated responses for the ratios At/t!^ = 1. The inset is a 
log-linear plot of averaged PDF of the scaled local integrated 
responses for the ratios At/tf^ — 1, the tail are well fitted by 
a Gaussian function ea;p(0.955 -I- 0.00881a; — 5.69613x'^). 



function collapse onto a master curve when plotted as 
tZC{t^ + Ai,U) and tZx{t^u + At,t^) versus At/t J. 
By minimizing the value of the difference bet-ween the 
master curve and the data we find that the best data 
collapse is obtained for the folio-wing aging exponents: 
a-f /3 = 0.34 ±0.05 and (3 = 0.48 ±0.05. We find (Fig. Ej) 
that the scaling functions satisfy the follo"wing asymptotic 
behavior: 



foiy) ^ fAiiy) y, 



(3) 



in agreement "with the fact that the motion of the parti- 
cles is diffusive at long times, (Aa;^) ~ At, and the exis- 
tence of a -well-defined mobility, respectively. Therefore 
at long times, both the correlation and response functions 
display the same po-wer la-w decay: 



- At,i^) 
At,t^) ' 



tJ^+^^At. 



(4) 



The result C(t^ + At^ty,) - + Ai,t^,) - i;^"-^^ 
indicates D{tyj) ^ M{t-u]) ~ tw^'^"^ ■ For short times, the 
MSD scaling function crosses over to a sub-diffusive be- 
havior of the particles. We obtain, {x^{tw + At.tw)) ~ 
t-"(At/tf,)0-3 = t-0-004AtO-3. Since the trapping time 
corresponds to the size of the cages denoted by q[tw), 
we can determine the cage dependence on t^ as q{tw) ^ 
t~" ™^. The resulting exponent is so small that we can 
say that the cages are not evolving with the waiting time, 
within experimental uncertainty. Furthermore, the scal- 
ing ansatz of Eq. |5J indicates that the relaxation time 
of the cages scales as T(t„) t^ since this is the time 
when the sub-diffusive behavior crosses-over to the long 
time diffusive regime. 

We find that the scaling forms of the correlations and 
responses are not consistent with p-spin models. Based 
on invariance properties under time reparametriza- 
tion, spin-glass models predict a general scaling form 
Cag(i, iu.) = Cag{h{t)/h{tw)) , where h{t) is a generic 
monotonic function We find that the scaling of 

our observables cannot be collapsed with the ratio 
h{t)/h{tu]). The scaling with h{t)/h{tu,) is expected for 
system in which the correlation function saturates at long 
times the other hand, our system is diffusive, 

and the studied correlation function is not bounded. In- 
deed, similar scaling as in our system has been found 
in the aging dynamics of anther unbounded system: an 
elastic manifold in a disordered media The sugges- 
tion is that this problem and the particle diffusing in a 
colloidal glass may belong to the same universality class. 
Furthermore our results can be interpreted in terms of 
the droplet picture of the aging of spin glass, where the 
growth of the dynamical heterogeneities control the ag- 
ing. 
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V. LOCAL FLUCTUATIONS OF 
AUTOCORRELATIONS AND RESPONSES. 

Previous work has revealed the existence of dynamical 
heterogeneities, associated with the cooperative motion 
of the particles, as a precursor to the glass transition as 
well as in the glassy state 0, 0, 0, ■ Instead of the 
average global quantities studied above, the existence of 
dynamical heterogeneities requires a microscopic insight 
into the structure of the glassy. Earlier studies focused 
mainly on probability distributions of the particles dis- 
placement near the glass transition. More recent analyt- 
ical work in spin glasses 16] shows that the probability 
distribution function (PDF) of the local correlation P{C) 
and the local integrated response P{x) could reveal es- 
sential features of the dynamical heterogeneities. 

Here we perform a systematic study of P{C) and P{x) 
in sample A, and the resulting PDFs are shown in Fig. |S1 
The scaling ansatz of Eq. (O implies that P{C) and P{x) 
should be collapsed by rescaling the time At by and 
the local fluctuations by f^. Indeed, this scaling ansatz 
provides the correct collapse of all the local fluctuations 
captured by the PDFs, as shown in Fig. \^ and O; 
for P{C) and in Fig. EJi for P(x). 

The PDF of the autocorrelation function displays a 
universal behavior following a modified power-law 



t-^P{C)^{tZC + Co)-\ (5) 

where Co and A only depend on the time ratio At/tf^. 
For the smaller values of C (C < t~"Co), the existence 
of a flat plateau in P(C) indicates that the tracers are 
confined in the cage. For larger values of C, the salient 
feature of the PDF is the very broad character of the 
distribution, with an asymptotic behavior P{C) ~ C^^. 
This large deviation from a Gaussian behavior is a clear 
indication of the heterogeneous character of the dynam- 
ics. Furthermore, the exponent A decreases from 2.6 to 
1.4 with the time ratio At/tf^ ranging from 10 to 60. We 
notice that A = 2 corresponds to the crossover between 
the short-time and long-time regime in Fig. 0] where 
At/tf^ w 40. The significance of A = 2 is seen in the in- 
tegral / P{C)CdC. For A > 2 {At/ti < 40) the plateau 
dominates over the power law tail in the integral and the 
dynamics is less heterogeneous. For A < 2 {At/t^ > 40) 
the power law tail dominates and this regime corresponds 
to the highly heterogeneous long-time regime. 

On the contrary, P{x), shown in Fig. displays 
a different behavior. The fluctuations are more narrow 
and the PDF can be approximated by a Gaussian. This is 
consistent with the fact that we did not find cage dynam- 
ics for the global response. Moreover, numerical simula- 
tions of spin-glass models |16|| seem to indicate a narrower 
distribution as found here. 



VI. SCALING BEHAVIOR OF LOCAL 
FLUCTUATIONS 

The local correlation function C and local integrated 
response function x studied here are calculated from each 
individual tracer trajectory. In a condensed colloidal 
sample, the tracer trajectories are always confined at 
a local position. Therefore, the correlation C and re- 
sponse X foi' ^n individual particle can be regarded as 
the coarse-grained local fluctuations of the observables as 
investigated in 16] . 

Furthermore, in order to improve the statistics of our 
results, the PDF P{C) is calculated not only for all the 
tracers, but also over a time interval much smaller than 
the age of the system. In practice, we calculate the i— th 
tracer's local correlation function as Ci(tu, -I- At,tw)- At 
a given At, we open a small time window [At — ts, At+tg] 
{tg ^ At) and count all the Ci{ty^ + Atj,t^) with 
Atj e [At -ts,At + ts] into the statistics of P(C). The 
calculated P{C) is a mixture of the local and the tempo- 
ral fluctuations of the observables. Similar technique is 
performed to calculate P{x)- 

Next we derive the scaling law for the PDF of the lo- 
cal correlation function. Let us first recall the scaling 
behavior of the global correlation function in Eq. lO, 



tZC{t^+At,t^) = foiAt/tJ), (6) 

where we add a bar to C to distinguish the global correla- 
tions from the local C. The average is taken over all the 
tracer particles. We can rewrite the global correlation 
function as the integration of the PDF P{C) as: 



C = J P{C)CdC. (7) 

Furthermore, we obtain the relation of fn and P{C) 
by substituting Eq. Q into Eq. ®: 

fD{At/tJ) ^tZC = J t-^P{C){tlC)d{tlC). (8) 

For a given At/tw^ , fo is equal to a constant, and Eq. 
© requires that t~°'P{C) should only depend on f^C. 
In other words, t~"P{C) is a function of At/tyj^ and 
t",C. Then we define Fd as: 

FD{tlC,At/ti)^t-'^PiC). (9) 

A similar formula can be obtained for P{x)- Eventually, 
we obtain the scaling ansatz of P{C) and P(x) shown in 
Fig.El 

P{C) =tlFD{tlC,At/ti), 

Fix) =t^FMitZx,At/ti), ^'^^ 
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where the universal functions F]j{x, y) and Fi\j(x, y) sat- 
isfy 

J FD{x,y)dx ^ foiy), ^-^x 
J FMix,y)dx = fMiy)- ' ' 



VII. STUDY OF THE PDF OF THE 
AUTOCORRELATION FUNCTION 

The PDF of the autocorrelation function follows a 
modified power law t^"P(C) oc (i",C + Cq)"'^, as we 
see in Figs. EId andEt. From Eq. we have 

tic^ r^^\-"p{c){tic)d{tic) 

_ ^^-\x + Co)-^xdx ^^^^ 
~ ^^-\x + Co)-^dx ■ 

The cutoff Ccut {Ccut ^ C'o) is introduced to make the 
integral converge and we always take A > 1, then 

cc=^[i-(A-i)(^+ir^]. (13) 

For A > 2, the last term in Eq. is negligible and 
t2,C « Co/(A — 2) mainly depends on the short-time pa- 
rameter Co. For A < 2, we have f^C « t5iC'o(%f )^"^ 
and the long-time parameter Ccut dominates. 

Following the previous discussion of P{C), we define 
the i-th tracer's local correlation function as: 

C,(i^ + At, tu,) = {Ax{tu, + At, tu,f),/2, (14) 

where the average is calculated for only the i-th 

tracer's trajectory by opening a small time window [At — 
ts,At-\- ts] [ts <IC At) at a given At, and counting all the 
Ci{tyj + Atj,tui) with Atj e [At - ts,At + ts] into the 
statistics of P{C). We should note that if there is no 
average P{C) can be reduced to a simple form: 

P{C) - P{{Ax^),/2) ^^^^ P{Ax^l2). (15) 



This form can be further reduced since 

P(AxV2)d(Aa;V2) = P{Ax)dx, (16) 
therefore at the limit of the average 0, 

P(C) = P{{Ax^),l2) ^•••^'""> P(AxV2) = P(Aa;)/Aa;. 

(17) 

Therefore our P{C) is related to the probability P(Aa;) 
usually studied in previous works ^TJ. The reader with 
enough endurance to have reached this point in the pa- 
per may realize the apparent contradiction between our 
power-law scaling and the typical stretched exponential 
behaviour found in other studies of local fiuctuations. In 
this regard we may first say that previous work consid- 
ered the supercooled regime while we work in the glassy 
regime. Thus, both regimes may show different type of 
scaling. When we analyze the scaling of P(Ax) in our 
glassy system we find that it can be approximately fit- 
ted by both a broad tail power law behavior for large 
Ax, consistent with equation [SJ but also by a stretched 
exponential, which is consistent with previous works 
for the supercooled regime. On the hand P{{Ax'^)i/2) 
can not be fitted by stretch exponential and it is only 
fitted by the power-law scaling. Therefore, we conclude 
that the power scaling should be the proper scaling for 
P(Aa;) and P(C) in the glassy regime. 



VIII. SUMMARY 

We have presented experimental results on the scaling 
behavior of the fluctuations in an aging colloidal glass. 
The power-law scaling found to describe the transport 
indicates the slow relaxation of the system. A univer- 
sal scaling form is found to describe all the observables. 
That is, not only the global averages, but also the lo- 
cal fluctuations. The scaling ansatz, however, cannot be 
described under present models of spin-glasses, but it is 
more akin to that observed in elastic manifolds in ran- 
dom environments suggesting that our system may share 
the same universality class. 
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